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Thermodynamics of a pseudospin-electron model without correlations is inves- 
tigated. The correlation functions, the mean values of pseudospin and particle 
number, as well as the thermodynamic potential are calculated. The calcula- 
tion is performed by a diagrammatic method in the mean field approximation. 
Single-particle Green functions are taken in the Hubbard-I approximation. The 
numerical research shows that an interaction between the electron and pseu- 
dospin subsystems leads in the p, = const regime to the possibility of the first 
order phase transition at the temperature change with the jump of the pseu- 
dospin mean value (S z ) and reconstruction of the electron spectrum. In the 
regime n = const, an instability with respect to phase separation in the elec- 
tron subsystem can take place for certain values of the model parameters. 
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1. Introduction 

The model considering an interaction of electrons with a local anharmonic mode 
of lattice vibrations has been used in the recent years in the theory of high-tempera- 
ture superconducting crystals. Particularly, such a property is characteristic of the 
vibrations of the so-called apex oxygen ions Orv along the c-axis direction of the 
layered compounds of the YBa 2 Cu 3 07-type structure [ [I], 0- An important role 
of the apex oxygen and its anharmonic vibrations in the phase transition into the 
superconducting state has already been mentioned [ f|, [5] and a possible connection 
between the superconductivity and lattice instability of the ferroelectric type in 
high-T c superconducting compounds has been discussed [ 0, [7J. In the case of a 
local double-well potential, the vibrational degrees of freedom can be presented by 
pseudospin variables. The Hamiltonian of the derived in this way pseudospin-electron 
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model has the following form [ |8| : 



(i.i) 



i ija 



and includes, besides the terms describing electron transfer (~ t^), the electron 
correlation ([/-term), interaction with the anharmonic mode (g-term), the energy 
of the tunnelling splitting (fi-term) and the energy of the anharmonic potential 
asymmetry (/i-term) in the single- site part 



Here, E gives the origin for energies of the electron states at the lattice site (E = 



In this paper, our aim is to obtain expressions for the correlation functions which 
determine dielectric susceptibility and the mean values of pseudospin and particle 
number operators, as well as the thermodynamic potential in the case of Q = and 
the absence of the Hubbard correlation U = 0. 

We perform the numerical calculations and investigate the mean values of pseu- 
dospin and particle number operators with the change of asymmetry parameter h 
(T=const) or temperature T (/i=const) for the cases of a fixed chemical potential 
value (regime /i=const) and a constant mean particle number (regime n=const). An 
analysis of the thermodynamic properties of the pseudospin-electron model in the 
case of absence of the electron correlation is also made. 

2. Hamiltonian and initial relations 

We shall write the Hamiltonian of the model and the operators which correspond 
to physical quantities in the second quantized form using operators of the electron 
creation (annihilation) at a site with a certain pseudospin orientation 



(1.2) 





Then, we obtain the following expression for the initial Hamiltonian: 



H 



hit 





+ 




where 



e = E + g/2, 



e = E -g/2 



(2.3) 
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are energies of the single-site states; H is a single-site (diagonal) term, H int is a 
hopping term. 

The introduced operators satisfy the following commutation rules: 
«jv} = 6ij5 aa >(- - S[), {af , a ja ,} = 0, 

a ja >} = SijS^i- + Sf), {at, a ja ,} = 0. (2.4) 

In order to calculate the pseudospin mean values we shall use the standard repre- 
sentation of the statistical operator in the form 

e -(SH = e -f>H z(p^ (2.5) 
a(P) = T T exp I f H int (r)dr } , (2.6) 



which gives the following expressions for (Si): 

(S?) = 7jrLr-(Sr*Wh = (Sfm)o- (2-7) 
Here, the operators are given in the interaction representation 

A(r) = e rHo Ae- rHo , (2.8) 

the averaging (. . .}o is performed over statistical distribution with the Hamiltonian 
H , and the symbol (. . .)§ denotes separation of connected diagrams. 

3. Perturbation theory for pseudospin mean values and a diagram 
technique 

Expansion of the exponent in ( |2.6|) in powers of H int ( [2.2|) leads, after substi- 
tution in equation fl2.7| ), to an expression that has the form of the sum of infinite 
series with terms containing the averages of the T-products of the electron creation 
(annihilation) operators fl2.1|) . The evaluation of such averages can be made using 
the Wick theorem. 

In our case this theorem has some differences from the standard formulation. 
Namely, each pairing of operators (|2.1|) contains operator factors, i.e. 



ai (T')at(T)= g(r' - r)5 lo P+, 5 l (r')5+(r)= g(r' - r)5 lo Pr, (3.1) 



a o( T ) a i( T ')= -9( T ' ~ r )hoP t + , 5+(r)a<(r')= -g(r' - r)^^ . 

Finally, this gives the possibility to express the result in terms of the products of 
nonperturbed Green functions 

9l ° { ] ~ (Ka+}) " e ^t-a + e^)- 1 ,r'>r, ^ 
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9io{T - r 



(T r a i (r)a+(r / )) 



( 1 + e -/3e)-l )T>r ' 

({^a+}) -(1+e^)- 1 ,r'>r. 

<7i (-r - r') = g{r - r')5 io , g io (r - t') = g(r - r')^, 



and averages of a certain number of the projection operators 



1 

2 



= - + si 



(3.3) 



Let us demonstrate this procedure for one of the terms which appear in the fourth 
order of the perturbation theory for {Sf)\ 



p p p p 
J dn J dr 2 J dr 3 /dr 4 ]TJ ^Aih^U 



P P P P 

«3J3 

(3.4) 

(I () o o *2J2«3J3 

x (TrSfaf (r 1 )a i (r 1 )S+ {r 2 )a h (r 2 )a+ (r 3 )h j2 (r 3 )a+ (r 4 )a j3 (r 4 )) . 

The stepwise pairing of a certain operator with the other ones gives the possibility 
to reduce expression (|3.4j) to the sum of the averages of a smaller number of operators 

(r T £jV( r i) a i( r i)2j (T 2 )a h (r 2 )a+ (r 3 )S j2 (r 3 )a+ (r 4 )a j3 (r 4 )) 



(T r 5f afiT^ajir^af {r 2 )a h (r 2 )a+ (r 3 )a i2 (r 3 )a+ (r 4 )a i3 (r 4 ) ) 



+ (T T 5f a+(ri)a J (ri)5+ {T 2 )a h (r 2 ) aj(r 3 )5 ja (r 3 )a+ (r 4 )a i3 (r 4 )) 
= -9ih{n ~ ^)(T r SfP+a i (ri)a+(r 2 )a il (r 2 )a+(r 3 )a i2 (r 3 )a+(r 4 ))o (3.5) 

- QijAn ~ ^)(T r 5fP+a j (ri)a+(r 2 )a+(r 3 )a j2 (r 3 )a+(r 4 )a i2 (r 3 )) . 



The successive application of the pairing procedure for ( |3.5|) leads, finally, to 

~9ih ( r i-r 2 )^ li2 (T2-T 3 )^ 3i (T4-T 1 )^ 1 (r 3 -r 2 ) (T T £f P/iJ.Pr j p+) 
-&ia ( r i-r 4 )^ li2 (r 2 -r 3 )^ 2i (r 3 -ri)^ 3il (r 4 -r 2 ) (T T S?P J + P J + 1 PrP+) 
+9ih(n-n)g hj2 (T2~T 3 )g i2jl (r 3 -T 2 )^ 

We introduce the diagrammatic notations 

1— D^O-l' " (^i' P i + + 9ivP{) 

and diagrams 



O; 



(3.6) 
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which correspond to expression 

Expansion of ( |3.6| ) in semi-invariants leads to multiplication of diagrams (semi- 
invariants are represented by ovals surrounding the corresponding vertices with di- 
agonal operators and contain the <5-symbol on site indexes). For example, 




We shall omit diagrams of types 2, 3, 5, i.e. the types including semi-invariants of a 
higher than the first order in the loop (this means that chain fragments form single- 
electron Green functions in the Hubbard-I approximation) and also connection of 
two loops by more than one semi-invariant (this approximation means that a self- 
consistent field is taken into account in the zero approximation). 

Let us proceed to the momentum-frequency representation in the expressions 
for the Green functions determined on a finite interval 0<r</3 when they can be 
expanded in the Fourier series with discrete frequencies 



g(r) 



-T, 

R 



1 



9{r) 



(3.7) 



1 



own) 



The characteristic feature of the already presented diagrams and the diagrams cor- 
responding to other orders of the perturbation theory is the presence of chain frag- 
ments. The simplest series of chain diagrams is 

+ ... , (3.8) 



2n + 1 



-7T. 




D=e- 



where 



g[uj r , 



(P + ) , (P~) 



(3.9) 



iuo n — e \u) n — e 

and corresponds to the Hubbard-I approximation for a single-electron Green func- 
tion. The expression 

1 



9' 



[UJn 



tk 



(3.10) 
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in the momentum-frequency representation corresponds to the sum of graphs ( |3.8| ). 
The poles of function Gk(u n ) determine the spectrum of the single-electron excita- 
tions 

e It n{tk) = l(2E + t k ) ± l -s]g 2 + tt k (S*)g + tl. (3.11) 

Behaviour of the electron bands as a function of the coupling constant is pre- 
sented in figure |l|. One can see that there always exists a gap in the spectrum. The 
widths of subbands depends on the mean value of the pseudospin and in the case 
of strong coupling (g 3> W) the subbands' halfwidth is equal to W{\ ± (S z )) (W is 
the halfwidth of the initial electron band). 




2.0 § 



Figure 1. Electron bands boundaries (W = 0.4, (S 2 



0.2). 



Let us now return to the problem of summation of the diagram series for the 
mean value (Sf) taking into account the above mentioned arguments. The diagram 
series has the form 




The analytical expressions for the loop has the following form 



(3.12) 




2 ^ tjil 

1^ „-!(,,, \ - 



N j-^ g- 1 ^) -t k \iuj n - e iu n - e 



p+ p: 
+ - 



(3.13) 
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where we used the notations 



2 tj 1__ 



2 tj 1 

a2 ~ N(3^(g-^ n )-t k ) {iwn-e)- 



Using decomposition into simple fractions and summation over frequency we obtain 

2 r i 2 -\ r 

a ±=Jj Yl tk [^i n ( £ i(^)) + s i w ( £ n(*fc)) » ft2 =^ tfc [^2w(ei(*i.)) + B 2 n(e u (t k )) 



where 



and n(£ 



^ = £i(t fc ) - e 

£i(*fe) — £n(*fc) 

^ = £i(tfc) - £ 

£i(^fe) — £n(tfc) 

is a Fermi distribution. 



Bo 



£n(*fc) — £i(^fc) ' 

£n(ifc) — £ 
£n(*fc) — £i(*fc) ' 



1 + lf e 

The equation for (Sf) can be presented in the form 
(Sf) = (St) -(Sf/3( ai P+ + a^Pr))l 



+ -ASt0\a l P+ + a 2 P l -) 2 Yo-... = (S t 



2! 



7?0 - j 9(aiP++a 2 Pr)jc_ 



(3.14) 



Let us introduce 



where 



rr _ \^ rrMF 



//f F = ^ + a 1 ^ + + o ;2J F»-. 
Then, the analytical equation for (Sf) can be expressed in the form 

Sp(Sfe-^ MF ) 

/Wmf) 

1 + e-P* 



(Sf) - (Sf)uF - 



= ^tanh y^(h + a 2 - oti) + In ^ • ( 3 - 15 ) 

The difference 0:2 — cci corresponds to an internal effective self-consistent field acting 
on the pseudospin 



a 2 — Oi\ 



£ — £ 



N ^ £ r (t fc ) - £„(tfc) 



w(£n(tfe)) - n(£ : (t fc )) 



(3.16) 
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4. The mean value of the particle number 



The diagram series for the mean value (rij) (using the perturbation theory, the 
Wick theorem and expansion in semi-invariants) can be presented in the form 



( n i) 




(4.1) 




where 



(P a ) 



iu n - e L 



+ ... 




P a =(P + ;P-), e a =(e;e) 



□ =fii 



and the last term appears due to the pairing of the electron creation (annihilation) 
operators with the operator of the particle number. 



An analytical expression for (4.1) can be obtained starting from formulae (|3[ 



(4.2) 



Using decomposition into simple fractions and summation over frequency we can 
present the mean value (rij) in the form: 

(n<) = ^ J2 [ n fe(^)) + "M**))] - 2(P+)n(e) - 2(p-)n(e). (4.3) 



5. Thermodynamic potential 

In order to calculate the thermodynamic potential let us introduce parameter 
A G [0, 1] in the initial Hamiltonian 



such that H — > H for A = and H — > H + H int for A = 1 . 
Hence 

Z x = Sp(e"^) = Sp(e"^a A (/3)) = Z (<7 A (/3)) C 

where 



a x {(3) = T r exp |-A J H- mt (r)dr^ 



(5.1) 
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and 



(5.2) 



AQ X = Q X -Q = ~ln(a x ((3)) . 

Here Qq is a thermodynamic potential calculated with the single-site (diagonal) part 
of the initial Hamiltonian. 
Therefore, 



AQ = dX 



dA 



(5.3) 



For value dfl x /dX, we can immediately write the diagram series in the next form: 



0- 



dX 




+ 



+ 




+ 



(5.4) 



where 



— Xt, and also 



+ h<£i 




Expression (|5.3| ) can be presented in the form (using the diagram series ( |5.4|) ): 



Att 



N/3 



V / Xtigl(uj n ) — V — — -dA 

i 



xt h ^4^ 



dA 



-dA. (5.5) 
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The first term in expression ( |5.5| ) may be written in a diagram form as 




(5.6) 



Series ( |5.6|) describes an electron gas whose energy spectrum is defined by the total 
pseudospin field. This series is in conformity with the so-called one-loop approxima- 
tion. 

The second term in expression ( f5.5| ) can be integrated to the following diagram 
series 




(5.7) 

and appears due to the presence of a pseudospin subsystem. 

Finally, the diagram series for (3AQ can be written as a sum of expressions ( |5.6|) 
and ( |5.7| ), and the corresponding analytical expression is the following: 



An 



Nf3 



(cosh |e)(cosh |e) 



1 [ (3 1 + e _/3e 

-- ln cosh + a 2 - at) + ln - - 



1, ,(p, , 1 + e"^ 

+— in cosh < —h + ln — 

(3 1 2 1 + e-l 36 



(S z ){a 2 -at) 



(5.1 



Here, decomposition in simple fractions and summation over frequency were done. 

Then, since the thermodynamic potential is a function of the argument (S z ), let 
us check the consistency of approximations made for (S z ), (n) and thermodynamic 
potential Q. For this purpose let us derive the mean values (S z ) and (n) from the 
expression for the grand thermodynamic potential 



d(-fi) N 



\ n ( £ M) + n(e a (t h ))] - 2(P + )n(e) - 2(p-)n(e), 



dHO 



1 f (3 1 + e _/3e 

2 tanh j ^(^ + "2 - «i) + In - - 
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We thus obtain 

Therefore, the calculation of the mean values of the pseudospin and particle num- 
ber operators as well as the thermodynamic potential is performed in the same 
approximation which corresponds to the mean field one. 



6. Pseudospin, electron and mixed correlators 

In this section our aim is to calculate the correlators 

K^{t-t') = (TSar)S z m (r')r, 
KZ{r-r') = (TS[(r)n m (r')) c , 
KZ{r-r>) = (Tn,(r)n m (r')) c , 



constructed of the operators given in the Heisenberg representation with an imagi- 
nary time argument. 

Let us present a diagram series for the correlation function (in the momentum- 
frequency representation) within a self-consistent scheme in the framework of the 
generalized random phase approximation (GRPA) (which was applied in [ |9], |T0| 
where the magnetic susceptibility of the ordinary Hubbard model and t — J model 
was considered). In our case (the absence of the Hubbard correlation) this approxi- 



mation is reduced, because the so-called ladder diagrams [ [TT| with antiparallel lines 
disappear. 

We would like to remind that we have omitted diagrams including semi-invari- 
ants of a higher than the first order in the loop and also connection of two loops by 
more than one semi-invariant. 



(s z s z ) q = ctrcr 




(6.1; 



where we define 




p« = (P+ j p-) ; a = (0,1); e a = (e,e) 



p a , o = s' 
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Here, the first term in equation ( |6.1| ) takes into account a direct influence of the 
internal effective self-consistent field on pseudospins and is given by 





3! 



o o 






+ . 



(6.2) 



Series (3.2) means a second-order semi-invariant renormalized due to the "single- 
tail" parts, and is thus calculated by H MF . 

The second term in equation (|6.1| ) describes an interaction between pseudospins 
which is mediated by electrons (the energy of the electron spectrum is defined by 
the total pseudospin field). 

We introduce the shortened notations 



n 



Q 





(6.3) 



Solution of equation ( |6.1| ) can be written in the analytical form 

(s z s z ) q = 



i + E(-i)^E^(| - (ST) 



(6.4) 



where 



3^ = 4j y %2tkti+ q r a (k,u n )r fi (k + q, 



N 



n,k 



T a (k,Lu r 



(6.5) 



(6.6) 



(iu n -e a ) (1 -t h g{u n ))' 

Decomposition of function r a (k,u n ) into simple fractions and the subsequent eval 
uation of the sum over frequency leads to the next expression: 



£(_l)o*/3 



n 



^2 



N ^ [e : (t fe ) - £n{tk)][£i{th+q) ~ £u(tk+q)} 

x f nle^tk)] - nle^tk+g)] n[e u (tk)} - n[e u (t k+q ) 
n[gifa)j - n[e n (tk +q )} _ n[e n (t k )} - n\ej^k+q)\ 



£l(tk) — £it(tk+q) 



(6.7) 
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we finally obtain an expression for 



After substitution ( |6.7|) in equation 
(S z S z ) q . 

This formula for the uniform case (q = 0) can be rewritten as 



(S z S z } q=0 =(l/4-(S z ) 2 ) 
x/l- 



(6.8) 



(e-e? 



_£i(tk) ~ Su(tk)] 3 



tl{e-e? 



{n[e I (t fe )] - n[e„(tfe)]} 



+ 



2iV ^ [ £l (i fc ) - e n (i fc )]2 cos h 2 £sM cosh 2 £sM / v 4 



Expression (|6.8| ) can be obtained from the derivative d(S z )/d(/3h). This means that 
the mean values of the pseudospin and pseudospin correlators are derived in the 
same approximation. 

For a mixed correlator the diagram series has the form 



(s z n) Q = Co_b) + Co~6 




(6.9) 



where 




(6.10) 



II 




(6.11) 



P°T a {k,u n ). 



Solution of equation ( |6.10| ) can be written in the analytical form 

I=2(n(e)-n{e))(S*S') q . 
Here we start from formula (|6.4j) and from the next relation: 

(S z n) MF -(S z )(n) 



2(n(e)-n(e)). 



(6.12) 



(6.13) 
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The second term in diagram series ( |6.9| ) can be presented as 
2 

n' 



£l(tk) ~ £n(tk) 

nje^tk)} - nle^tk+g)] n^fa)] - n[e u (tk +q )} 

£l(tk) — £i(tk+q) £ l(tk) — £ ll[tk+q) 

n[e n (t k )) - n[ei(tjn-,)] n[e u (t h )) - n[e u (t h+q ))~ 



Let us introduce the shortened notations for expression ( 6.14j) 

II = (S'S*) q [®] q . 

In this way we obtain 



(S z n) t 



me) — me) 



{S*S") q +(S K S*) q [®]< 



(6.14) 



(6.15) 



(6.16) 



From our diagram series we can see: correlators containing the pseudospin variable 
S z are different from zero only in a static case. This is due to the fact that operator 
S z commutes with the Hamiltonian, being an integral of motion. 
For an electron correlator our diagram series has the form: 




(6.17) 



and only the last term is not equal to zero for non-zero frequencies. Let us consider 
series ( |6.17|) term by term. 

The first term in series ( 6.17|) may be written as 




(6.18) 



a, p 

After simple transformation we can obtain the next relation: 

1 / (P + ) IP-) \ {(nS z ) MF -(n)(S z )) 2 



(nn) mf - (n) - 



(P + ) | (P~) 
2 \ cosh 2 f cosh 2 f 



(P + )(P-) 



(6.19) 



This relation makes it possible to write immediately a simple analytical expression 
for series fl6.18| ) 



[2(n(e)-n(e))] 2 (S z S Z ) 



(P + ) 



(P-) 



q 2 I cosh 2 ^ ' cosh 2 ^ 



5{u). (6.20) 
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Analytical expressions for the //-term can be obtained starting from formulae 

(6.21) 



II={2[n(e)-n(e)){S*S z )M q \6(u) 
Using expression ( |6.16|) we can unite ( |6.21| ) and (|6.20|) 



/ + // 



2(»(,)-n(.)) (S .»>, + i M^W + ^l 



The diagram series for the fourth term in ( |6.17|) has the form 




6(u). (6.22) 



and can be written as 



IV 




Using formula ( 6.1 6| ) once more we unite the ///-term and the IV-term 



III + IV = (nS z ) q [@] q 5{uj). 
The last term can be presented in the form 





a , p 



(6.23) 



(6.24) 



(6.25) 



Let us write down the final formula for an electron correlator for the uniform (q = 0) 
and static (u = 0) case 



(nn) = 2(n(e)-n(e))(S z n) q=0 + 



(P + ) + (P~) 



2 \ cosh 2 f cosh 2 f 



(6.26) 



this - e) 



-E 



+ 



2N ^ £l (t h ) - e a (t h ) { cosh 2 cosh' , 

1 



2 gglftfe) 



(S z n) q=0 



2N V I cosh 2 £sM cosh 2 Mil 



1 J 1 1 

2 jrff + cosh 2 f 



The same result can be obtained from the derivative d(n)/(d/3/x). Thus, all our 
quantities: the mean values of the pseudospin and particle number operators, the 
thermodynamic potential as well as correlation functions are derived within the 
framework of one approximation which corresponds to the mean field approximation. 
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7. Numerical research in the \± = const regime 



In the investigation of equilibrium conditions we shall separate two regimes: 
/i=const and n=const. For the first regime the equilibrium is defined by the minimum 

= 0. 



T,/j,,h 



of the thermodynamic potential: [ jSgz\ 

The /x=const regime corresponds to the case when charge redistribution between 
the conducting sheets CuC>2 and other structural elements (charge reservoir, e.g., 
nonstoichiometric in oxygen CuO chains in YBaCuO-type structures) is allowed. 

The dependencies of the order parameter (S z ) on field h and temperature T at 
the constant value of the chemical potential are determined by equation ( 3.15|) . All 
the integrals in ( |3.15D can be calculated analytically at zero temperature (below, 
all the calculations will be performed for the rectangular density of states, but we 
would like to note that the similar behaviour can be obtained in the case of the 
semi-elliptic density of states). 

We shall present all our results for the case of zero temperature as well as for 
the case of non-zero temperature. 



a 1.5 




-0.5 0.0 



0.5 l-l 1.0 



-1.0 -0.5 0.0 



0.5 |^ 1.0 



Figure 2. Phase diagram ^ - h. Dotted and thin solid lines surround regions with 
S z = respectively. Thick solid line indicate the first order phase transition 
points, a) the case of zero temperature; b) T = 0.002 (W = 0.2; g = 1) . 

The phase diagrams which indicate stability regions for states with 
(S z ) = ±| are shown in figure for g ^> W. 

One can see two regions of the /i and h values where the states with (S z ) = ~ and 
(S z ) = — \ are both stable. In the vicinity of these regions the phase transitions of 
the first order with the change of the longitudinal field h and/or chemical potential 
\x take place and they are shown by thick lines on phase diagrams in figure ^. 

The field dependencies of (S z ) and Q near the phase transition point are pre- 
sented in figures |3], ^. Their behaviour in the cases of T = and T / with the 
change of the chemical potential is similar: S-like for the mean value of the pseu- 
dospin and a fish tail form for the thermodynamic potential. 
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Figure 3. Field dependency of (S z ) (W = 0.2, \i = —0.4, g = 1.0) for /i=const 
regime; T = 0.01 and T = 0. 



In the /i=const regime, the chemical potential can appear in the electron bands 
or out of them with the change of field h (dashed line in figure [5]), and in the vicinity 
of the phase transition point this results in a rapid change of electron concentration 
(dotted lines in figure §) due to a charge transfer from/to the reservoir (CuO planes). 
The widths of the electron subbands depend on the mean value of the pseudospin 
which results in the presented above behaviours. 

The phase transition point is presented by a crossing point on the dependence 
Q(h) (figure |j). At the same time this value is determined according to the Maxwell 
rule from the plot of function S z (h). 
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Figure 4. Field dependency of thermodynamic potential (W = 0.2, fx = —0.4, 
g = i.o); T = 0.01 and T = 0. 

With the temperature increase the region of phase coexistence narrows, and the 
corresponding phase diagram T c -h is shown in figure 0a. The tilt of the coexistence 
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Figure 5. Field dependency of electron bands boundaries (W = 0.2, \i = —0.4, 
g = 1.0); T = 0.01 and T = 0. 
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Figure 6. Field dependency of electron concentration (W = 0.2, fi 
g = 1.0); T = 0.01 and T = 0. 
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curve testifies to the possibility of the first order phase transition at a change of 
temperature with a jump of the pseudospin mean value. (The phase diagram T c -\i 
has a similar form). The existence of the shifted and tilted coexistence curve as the 
result of the local pseudospin-electron interaction was obtained for the first time in 



Y2\ for a pseudospin-electron model with a direct interaction between pseudospins. 



The analysis of the thermodynamic potential behaviour with the temperature 
increase (figure |7|b) shows that the lowest value of f2(T) corresponds to the jump of 
the mean value of the pseudospin (dotted lines in figure |8|) from the branch which 
corresponds to the low temperature phase to that of the hight temperature phase. 
The analysis of the (S Z S Z ) behaviour with the temperature decrease shows that 
the high temperature phase is stable up to zero temperature. This means that the 
vertical line on the T c -h phase diagram only once crosses the boundary of the phase 
stability. 

In figures the case when the chemical potential is placed in the lower subband 
is presented. There is no specific behaviour when the chemical potential is placed 
out of the bands. And if the chemical potential is placed in the upper subband our 
results transform according to the internal symmetry of the Hamiltonian: 



-H, h^2g-h, 



n 



2-n, 



S z 



-S z 
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Figure 7. a) phase diagram T c -h (g = 1, fi = —0.4); b) temperature dependence 
of the thermodynamic potential (W = 0.2, h = 0.22, g, = —0.4, g = 1). 



8. Numerical research in the n = const regime 

In the regime of a fixed value of electron concentration the first order phase 
transition with a jump of the pseudospin mean value accompanied by a change of 
electron concentration transforms into a phase separation. 

The dependence of the mean value of the particle number (or electron concentra- 
tion) on the chemical potential is one of the factors determining thermodynamically 
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stable states of the system. One can see the regions with d/x/dn ^ where states 
with a homogenous distribution of particles are unstable, which corresponds to the 
phase separation into the regions with different electron concentrations and pseu- 
dospin mean values (figures |9] and [Uj) . 

In the n=const regime the equilibrium condition is determined by the minimum 
of free energy F = Q + fiN. In the phase separated region the free energy as a 
function of n deflects up (figure |H]) and concentrations of the separated phases are 
determined by the tangent line touch points (these points are also the points of 
binodal lines which are determined according to the Maxwell rule from the function 
fi(n), see figure ^). 

The resulting phase diagram T-n is shown in figure |ll]. 
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Figure 8. Dependence of the mean value of the pseudospin and the pseudospin- 
pseudospin correlation function on temperature (W = 0.2, h = 0.22, fx = — 0.4, 

9 = 1). 



9. Conclusions 

Investigation of a pseudospin-electron model in the case of electron correlation 
absence was performed in the mean field approximation using the Hubbard-I ap- 
proximation for the calculation of a single-particle Green function. We presented the 
analytical consideration of our model and all the quantities were obtained within 
the framework of one self-consistent approximation. 

As the result of numerical investigations we have obtained: 

1) there is always a gap in the electron spectrum (with a change of the mean 
value of the pseudospin a reconstruction of the electron spectrum takes place); 

2) the possibility of the first order phase transition with a change of the longi- 
tudinal field h (as a consequence of this the S-like behaviour of the mean value of 
the pseudospin with a jump in the phase transition point (which corresponds to the 
inflected point on the dependence Q(h)) is obtained and at this point the concen- 
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Figure 9. Dependence of the chemical potential fi on the electron concentration 
n and the pseudospin mean value (S z ) for different h values (g = 1, W = 0.2, 
T = 0.01). 




Figure 10. Dependence of the chemical potential /i on the electron concentration 
n and deviation of the free energy from linear dependence for different T values 
( g = l,W = 0.2, h = 0.2). 
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Figure 11. Phase diagram T-n (g = 1, W = 0.2 h = 0.2) for a phase separated 
state. Solid line indicates binodal points, dashed line indicates spinodal points. 



tration rapidly redistributes between the conducting sheets CuC>2 and the charge 
reservoir (CuO planes) in YBaCuO-type structures); 

3) the phase coexistence curve is tilted from the vertical line, therefore, there 
exists a possibility of the first order phase transition with the temperature change; 

4) the high temperature phase is stable in the whole region of temperatures 
(figure |); 

5) in the regime n=const we have the regions with d^i/dn < 0, which corre- 
sponds to phase separation with the appearance of regions with different electron 
concentrations and different orientations of pseudospins. 

Analytical expressions for the mean values, thermodynamic functions and sus- 
ceptibilities of such a simplified pseudospin-electron model (U — 0, ft — 0) were 
obtained for the uniform phase, when (S?) = (S z ), and only the possibilities of 
phase transitions with uniform changes (q = 0) were analysed numerically. On the 
other hand, it is known that for certain parameter values the charge ordered phase 
can exist in a strong coupling limit of the pseudospin-electron model (U — > oo) [ 



13| which calls for the consideration of the possible superstructure orderings in the 



opposite limit of U = and will be the subject of our further investigations. 
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